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Integrate and fire oscillators are widely used to model the gener- 
ation of action potentials in neurons. In this paper, we discuss small 
noise asymptotic results for a class of stochastic integrate and fire 
oscillators (SIFs) in which the buildup of membrane potential in the 
neuron is governed by a Gaussian diffusion process. To analyze this 
model, we study the asymptotic behavior of the spectrum of the fir- 
ing phase transition operator. We begin by proving strong versions 
of a law of large numbers and central limit theorem for the first 
passage-time of the underlying diffusion process across a general time 
dependent boundary. Using these results, we obtain asymptotic ap- 
proximations of the transition operator's eigenvalues. We also discuss 
connections between our results and earlier numerical investigations 
of SIFs. 



1. Introduction. The integrate and fire oscillator is widely used to 
model the behavior of the membrane potential in a neuron. Since its intro- 
duction by Lapicque [6] in 1907 it has been studied by many authors, see in 
particular Stein [13] and Knight [9]. The neurobiological derivation of the 
model is described in Tuckwell [18, 19] and a recent review of activity in the 
field appears in Burkitt [2, 3]. 

In this paper we use the following model for the stochastic integrate and 
fire oscillator (SIF). After starting at some time to the membrane potential 
evolves according to the stochastic differential equation 

(1) dXf = {-iXl + I{t))dt + e dWt 

until it reaches a (time-dependent) threshold level g{t) at time 

(2) rf = mf{t > to : Xf = g{t)}. 
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At the hitting time rf the membrane potential discharges, producing a volt- 
age spike, and resets at a lower value 

(3) Xf^.^^^=h{t). 

For t > (rf)"*" the process follows the SDE (1) until the second hitting 
time t| = inf{t > rf : = g{t)}, and so on, yielding a sequence {r^ : n > 
1} of hitting times. 

Here Wt is a standard one-dimensional Wiener process, and e > de- 
termines the intensity of the noise in the integrate and fire oscillator. The 
input function I{t), the threshold function g(t) and the reset function h{t) 
are deterministic functions and will be regarded as given as part of the 
problem. The parameter 7 > gives the rate of leakage of current across 
the membrane. The terms "leaky" and "non-leaky" are sometimes used to 
describe the cases 7 > and 7 = respectively. The issue is to provide a 
concise description of the distribution of the random sequence {r^ : n > 1}. 
It is of particular interest to describe how the distribution of the sequence 
{r^ : n > 1} responds to changes in one or more of the functions /, g and h. 

If all three functions are constant, then the inter-spike intervals t^+i — 
form an independent, identically distributed sequence of random variables. 
A more interesting situation occurs when one of the functions undergoes a 
periodic modulation. In many applications the input function is taken to be 
of the form I{t) = Iq + Iisinuit. In other cases the threshold is taken to 
be of the form g(t) = go + gi sinuit. In this paper we will make the general 
assumption that the three functions / and g and h all have the same period. 
Without loss of generality we will assume that the period is 1. Then the 
sequence of firing phases 

©n = < mod 1 

determines a Markov chain {G^ : n > 1} on the circle S = R/Z. 

When e = the process X^ is given by an ordinary diff'erential equation. 
Therefore the hitting times are given by r^+i = /(t^) for some deterministic 
function / satisfying /(t + 1) = /(i) + 1, and the firing phases are given by 
Qn+i ~ /(©n) where f{6) = f{9) mod 1. In settings where either the input 
function I{t) or the threshold function g{t) is of the form A + B sin27rt (and 
the other two functions are constant), the dynamical system on S generated 
by iterating / has been studied by Rescigno, Stein, Purple and Poppele [11], 
Knight [9], Glass and Mackey [5] and Keener, Hoppensteadt and Rinzel [8]. 
Of particular interest are the regions in the {A, B) parameter space giving 
rise to phase-locked behavior, and the bifurcation scenario as A and B are 
varied. 
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When £ > the deterministic hitting time function / is replaced by the 
first passage-time density function 



and / is replaced the projection p{9\6q), say, of p'^(t|io) onto the circle S. The 
behavior of the Markov chain {0^ : ?^ > 1} may be studied via its transition 
operator given by 



for (j) in the class B{S) of bounded measurable functions on S. For any 
e > the transition densities p^{0\9q) are bounded away from zero, so that 
the Markov chain {0^ : n > 1} is uniformly ergodic and has a unique 
stationary probability distribution. The compact operator captures the 
essential dynamics of 0^, and hence its spectrum is of primary interest in 
quantifying the transient and asymptotic behavior of the system. 

In a sequence of papers Tateno [14, 15] and Tateno and Jimbo [16] consider 
the effect of small noise on the deterministic bifurcation scenarios considered 
earlier. The papers [14, 15, 16] contain numerical calculations of the leading 
eigenvalues of the transition operator . These calculations suggest a qual- 
itative change in the small noise behavior of the leading eigenvalues near the 
location of the deterministic bifurcation. The calculations in [14, 15, 16] in- 
volve numerical approximations in two places. Firstly, since there no explicit 
formula for the first-passage density p{t\to) except in a few special cases, nu- 
merical techniques are used to solve an integral equation for p{t\tQ), following 
the method proposed by Buonocore, Nobile and Ricciardi [1]. Secondly, the 
circle S is replaced by a finite set of points. Thus the operator acting on 
B(S) is approximated by a finite-dimensional stochastic matrix. 

In this paper we obtain rigorous results on the asymptotic behavior of 
the spectrum of the operator T'^ as e ^ 0. The first main result gives a 
Gaussian approximation for the first-passage density p^{t\tQ) as £ 0. This 
result does not use the assumption of periodicity, and is valid for any 
functions I(t) and g{t) and starting position X^{tQ) = xq < g(to), under the 
condition that the deterministic hitting time is finite and that the determin- 
istic trajectory crosses the threshold transversally. For details see Section 2 
and especially Theorem 1. This result can be applied in the periodic setting 
to show that the Markov chain {G^ : n > 1} can be well approximated by 
small Gaussian perturbations away from the deterministic mapping /. The 
estimate in Theorem 1 is sufficiently strong that the techniques in Mayberry 



p'{t\to) := 7^F« < t|r; 



i=to) 
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[10], which deals with small Gaussian perturbations of circle maps, can be 
applied here also. In the simplest case where the deterministic mapping / 
is continuous and has one stable fixed point Og attracting all orbits except 
the one started at one unstable fixed point 9u, the limiting eigenvalues of 
can be calculated explicitly in terms of f'{Os) and f'{6u), see Theorem 2 
in Section 3. This result can be extended to the case where / is continuous 
and phase- locked, see Remark 3. In many examples of SIF the determin- 
istic mapping / has a finite set of discontinuities, and this case is treated 
in Section 4. The main result. Theorem 3, deals with the case where / is 
phase-locked and where the discontinuities are well away from the phase 
locked orbit. (This rather vague assertion is made precise in condition (D3) 
of Theorem 3.) 

Tateno and Jimbo [16] consider the leaky SIF with constant input I, 
periodically modulated threshold g{t) = 1 + A;sin27r and constant reset level 
0. Using a 100 x 100 stochastic matrix in place of the operator T^, they 
produce plots of the leading eigenvalues for various small values of the noise 
intensity e. In Section 5, we indicate how, in the phase-locked setting, our 
results may be applied to give a theoretical interpretation of the e — > 
behavior seen in some of the figures of [16]. 

Finally, Sections 6, 7 and 8 contains the proofs for the results in Sections 
2, 3 and 4 respectively. 

2. First passage times. The results in this section do not use peri- 
odicity, and are valid for general functions I{t) and g{t), and for all 
7 > 0. 

Denote by P*0'^o ^j^g law of the diffusion process {X^ : t > Iq} satisfying 
(4) dXf = (-7Xf + I{t))dt + e dWt 

with initial condition X^^ = xq. For xq < g(io) define the first passage time 

= inf{t > to : = g{t)} 

of the process Xf across the threshold g{t). In this section we will consider 
the behavior of the distribution of , and in particular its density function 

p'{t\to,xo) = ^F'^'^"{T' <t), 

as e — > 0. 

Let (^{t) = £^{t\tQ,xo) denote for t > to the solution to the noise free 
(e = 0) equation (4) with initial condition ^(to) = < g{to). Thus 

e(t|to,^o) =e-^(*-*«)xo+ re-T(*-^)/(s)ds. 

J to 
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For xo < 9(^0) let 

/(to, xq) = mf{t > to : C{t\to, xq) = g{t)} 

denote the deterministic hitting time. If f{to,XQ) < 00, define 

m(to, xq) = -jgifito, Xo)) + /(/(to, xq)) - g'{f{xo, to)). 

Thus m(to,xo) measures the difference in slopes when the deterministic so- 
lution ,^(t|to,xo) first meets the threshold g{t). Since xo < g{to), the deter- 
ministic solution hits from below and so m(to,xo) > 0. The deterministic 
solution crosses the threshold transversally if and only if m(to, xo) > 0. Now 
define the set 



g = {(to,xo) : Xo < g{to) and /(to,xo) < 00 and m(to,xo) > 0}. 

of initial conditions (to,xo) for which the deterministic trajectory crosses 
the threshold function transversally at the finite time /(to,xo). 

Proposition 1. Q is an open set and f G C^(^,M). 

The proofs of the results in this section can be found in Section 6. Our 
next result gives a uniform bound on the deviation of away from the 
deterministic crossing time /(to,xo). 

Proposition 2. Let G be a compact subset of Q. Then for any 5 > 
there are constants Ms and Ks such that 

P*o,xo(|^. - /(to,xo)| >6)< M,ee-^^/^' 
for all (to, Xo) G G. 



For (to,Xo) G Q define 

(1 - e-27(/{to,xo)-to))/2^ if 7 > 0, 



a (to,xo) 



/(to,xo)-to if 7 = 0. 



Notice that e'^a'^{tQ, xq) is the variance of at the moment t = /(to,xo) of 
noise-free intersection. Define 

/(-N 2/+ \ cr^(to,Xo) 
(5) C7^(to,Xo) = r 

m^(to,xo) 
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and 



so that PtI^I^Oi^jo) is the density at t of a N{0,a'^{tQ,xo)) normal random 
variable. With this notation in hand, we can state the main result of this 
section. 

Theorem 1. Let G be a compact subset of Q . Then there exist finite 
positive constants 5, ai, K and Eq (depending on G) so that 

(6) sup |ep"(/(to, xq) + u|to, xo) - Pr(n/e|io, xo)\ < ifee-^'/^^'^i 
for all £ < Eq, \u\ < 6. 




Fig 1. First passage-time pdf for dXt — {—Xt + 2)dt + .IdWt, Xo = .5, across g{t) — 
1. Dashed line: Gaussian approximation of Theorem 1 with mean ln(1.5) ~ .4055 and 
standard deviation -^5/1800 ~ 0.0527. Solid: Numerical approximation obtained by solving 
integral eguation from [1] 
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Remark 1. Under P*0'^o the centered and scaled hitting time (r"^ — 
f{tQ,xo))/e has density ep^{f{tQ,XQ) + et\tQ,xo) at t. Putting u = et in 
Theorem 1 gives the result that under P*0'^o 

(7) ^^-^^^^^N{0,aUto,xo)) as 6^0, 

or more informally 

N{f{to,xo),e'^a^{to,xo)) as e ^ 0. 

Figure 1 shows the densities of and N{f{to,xo),e'^ar{to,xo)) for an ex- 
ample with e = 0.1. Of course the result in Theorem 1 is much stronger 
than (7) since it gives locally uniform convergence of densities, rather than 
just convergence in distribution. This extra strength will be important for 
the spectral analysis results in the next two sections. 

Remark 2. In equation (5) the difference in slopes m{to, xq) at the point 
of deterministic intersection is used to convert the variance cj(to,a;o) in the 
spatial dimension into the variance ar{to,XQ) in the temporal dimension. A 
heuristic observation of this conversion factor appears in Stein [13]. 

3. Transition Operator for the SIF. We now return to the setting 
of SIFs. The functions I{t), g{t) and h(t) are assumed to be and periodic 
with period 1. Also h{t) < g{t) for all t, and 7 > 0. The deterministic 
solution is 

e(t|to, xo) = e-^(*-*«)xo + f e-^(*~^)/(s) ds. 

Jto 

If 7 = 0, then the assumption I{s)ds > implies that S,{t\to,xo) 00 
as t ^ 00 and so the deterministic hitting /(to,a^o) is finite. If 7 > then 
^(t\to,xo) — (,{t) — > as t — > 00 where 



1 - e-T Jt~i e-z-lJo 

In order to ensure that crossings for the noise free system occur in a nice 
enough fashion we impose the following conditions on the functions I{t) and 



max{^(t) - g{t) :0<t<l}>0 if7>0 

/^/(s)(is>0 if 7 = 

Jo 



(A) 

(B) --fg{t) + I{t) - g'{t) > for ah t 
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Note that (A) implies that f{tQ,XQ) < oo whenever xq < g{to), and then 
(B) imphes that m{tQ,XQ) > 0. Therefore Q = {(tQ,XQ) : xq < g^to)}, and 
we can apply our results on first passage densities from Section 2 using the 
compact set G = {{s,h{s)) : < s < 1}. We always take xq = h{to), so we 
write p^(t|to) = p'{t\to,h{to)), /(to) = f{to,h{to)), a^ito) = a^^{to,h{to)), 
etc. Proposition 1 implies that / G C^(M), and clearly f{t + 1) = f{t) + 1 
for all t G M. Theorem 1 implies that when e > is small, we have the 
approximation 

(8) < ~ + e(^rK-l)Xn 

where {xn '■ n- > 1} is & sequence of independent A^(0, 1) random variables. 

The sequence of firing phases {0^ : n > 1} is a Markov chain on the circle 
S = M/Z with transition density function 

(9) fie\9o)=J2p'(^ + ^\^o) 

for all 9, 00 G S. Looking at (8), we may expect that the study of {6^ : n > 1} 
should be similar to the study of the chain 

(10) = f{Y^_,) + eaAY^_,)xn mod 1. 

Spectral properties of the transition operator for Markov chains of the form 
(10) were developed in Mayberry [10]. Here, we prove a similar result for 
the transition operator of the chain {0^ : n > 1}. 

Theorem 2. Let {G^ : n > 1} be the sequence of firing phases for a 
period 1 SIF with input, threshold and reset functions I{t), g{t) and 
h{t). Assume the conditions (A) and (B). Suppose that the deterministic 
phase return map 

f = f mod 1 

has a stable fixed point 6s and an unstable fixed point 9^, and that f'^{9) — > 9s 
for all 9 € '§i \ {9u}. Let denote the transition operator for {Q^ ■ 1}. 
Then for any r > and e sufficiently small, we can write = T^^ + T^p 
where \\T^p\\oo < r and any eigenvalue of T^p with modulus greater than r 
is of one of the two forms c" + 0{e) or |cu|~"'^c~"' + 0{e) for some n > 0, 
where Cg = f'{9s) and Cu = f'{9u). 

The proof of this result can be found in Section 7. The result says that 
and |cm|~^c~" are the limiting eigenvalues of in the sense that for 
any r > 0, has sequences of r-pseudoeigenvalues which converge to c" 
and \cu\~^c~'"' as e — > (see Trefethen and Embree [17] for definitions of 
pseudoeigenvalues) . 
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Remark 3. Theorem 2 can be extended to the case where / has pe- 
riodic orbits. Consider an orbit {61,62, ... ,0^^} of period k > 1 and let 
c = /'(^i)/'(^2) • • • fiGn) denote the product of the derivatives of / along 
the periodic orbit. If the orbit is stable, so that |c| < 1, then it contributes 
limiting eigenvalues (c")"^/'^ for n > 0. If the orbit is unstable, so that |c| > 1, 
then it contributes limiting eigenvalues (|c|~^c~'^)^/'^. Here all the Kth roots 
are included as limiting eigenvalues. The proof in this more general setting 
combines the method of proof of Theorem 2 with techniques used in the 
proof of [10, Theorem 1], and is left to the reader. Moreover, the methods 
used in [10, Theorem 2] can be applied here to give information about the 
associated eigenfunctions. 

4. Discontinuous Case. Tateno and Jimbo [16] consider several cases 
of SIFs in which / is well defined, but Assumption (B) fails and / is 
discontinuous at some 6* G §. In this section, we will discuss extensions 
of Theorem 2 to this situation. As before, for /i(to) < g{to) we define 
/(to) = inf{t > to : iiA'^o) = g{t)} < but now we also define /*(to) = 
inf{t > to : '^(t|to) > g{t)}. Thus /(to) is the time of first hitting of the 
threshold, and /*(to) is the time of first crossing of the threshold. We keep 
(A) unchanged, but replace condition (B) with 

(B') There is a finite set D C [0, 1) (possibly empty) such that — 75(/(to)) + 

lifito)) - g'ifito)) > for all to G [0, 1) \ D. 
(C) For each to £ D either 

(i) f*{to) = f{to); or else 

(ii) /*(to) > /(to) and git) > e(t|to) for /(to) < t < /*(to). 

Notice that (A) implies that /(to) < /*(to) < 00 for all to, and that (B') 
implies that Proposition 2 can be applied to any compact subset of [0, 1] \ D. 
In case (C')(i) the deterministic trajectory ^(t|to) crosses the threshold g{t) 
at t = /(to), and / is continuous (but not differentiable) at to; and in case 
(C')(ii) ^(t|to) touches g{t) at t = /(to) and then does not intersect again 
until it crosses at time t = /*(to), and / is discontinuous at to. Examples of 
the behavior described in (C) are given in Figure 2. 

We begin with the following extension of Proposition 2. Proofs for the 
results in this section can be found in Section 8. 

Proposition 3. Suppose that f satisfies (A), (B') and (C), and D / 
0. Then for all 6 > there exist 6 > and K, M such that 

P*(d(r^{/(to),r(to)}) >5)< sMe-""'' 
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whenever \t — IqI < 5 for some to & D. 

The assumptions in our final theorem are not the most general ones possi- 
ble, but the result is sufficient to treat the examples in the next section and 
to give the reader the indication of how Theorem 2 (see also Remark 3) can 
be extended to discontinuous settings. The conditions may seem awkward, 
but we will see in the next section that they are easy to verify numerically 
in examples of interest. 

Theorem 3. Suppose that (A), (B'), and (C) are satisfied and in 
addition that f = f mod 1 satisfies the conditions 

(Dl) / has a periodic orbit P = {6i, . . . , 0^} in§ \ D for some k> 1 and 

\f'{e,)f'ie2)---f'ie,)\<i. 

(D2) r{e) ^ P as oo for all e £ S. 

(D3) /"^(-D) = for some £ > 1 and r{E) n D = ^ for < i < £ - 1, 
where E = f{D)Uf*{D). 

Let denote the transition operator for Q^. Then for any r > and e 
sufficiently small, we can write = Tf^ ~\~T^p where ||T;p||oo < CLnd any 
eigenvalue ofT^p with modulus greater than r is of the form (c")^/'^+0(e) for 
n > 0, where c = /'(0i)/'(02) • • • /'(^«) = ih'{Oi) for any i £ {1, . . . , k} . 

5. Examples from Tateno and Jimbo. We now apply our spectral 
results to the examples considered in [16]. There, the authors consider leaky 
SIFs with constant input I{t) = I, threshold g{t) = 1 + ksinlirt, and reset 
h(t) = 0. Since ^(t) = I/j the condition (A) becomes 

(11) I/-f>l-k. 
The transversality condition (B) is now 

-7(1 + k sin 2-Kt) + I - 2iTk cos 2-Kt > for all t 
or equivalently 

(12) 7/7 - 1 > fcy^47r2/72 + 1. 

If condition (B) fails, then there is a point of tangency on the threshold curve 
which leads to a point with the property described in (C')(i) or (C')(ii). 
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In any case, assuming (A) holds, we obtain /(to) = inf{t > to '■ h{t) = to} 
where 



If h'{t) > for ah t then / = h^^ is a smooth function. However, strict local 
maxima of h give rise to discontinuities in /. Moreover periodic orbits of the 
induced mapping / on S can be found using the facts that if h'^{to) = Iq mod 
1 then /^(to) = to mod 1 and (/'')'(to) = l/{h^y{to). The plots of / and 
the numbers in the examples below are obtained using the explicit formula 
for h. 

Tateno and Jimbo take 7 = 1/12.8 and various values of / and K. 

Example 1. Taking 1=1 and k = 0.1 gives the map / in Figure 3, 
which satisfies the conditions of Theorem 2. There is a stable fixed point at 
0.5622 with /'(0.5622) = 0.6142 and an unstable fixed point at 0.9379 with 
/'(0.9379) = 2.6898. All orbits starting at 9o ^ 0.9379 converge to 0.5622 so 
by Theorem 2 the limiting eigenvalues are {0.6142"' : n > 0} U {2.6898~"'~^ : 
n > 0} = {1, 0.6142, 0.3772, 0.3718, 0.2317, . . .} 

Example 2. Keeping 1 = 1 and increasing k to k = 0.35 gives the map 
/ in Figure 4. There is a discontinuity at 0.1178 with /(0.1178) = 0.8208 
and /*(0.1178) = 0.3946 so in the notation of Theorem 3, we have D = 
{0.1178} and E = {0.8208,0.3946}. Clearly /"^(D) = and E n D = Hi so 
that (D3) is satisfied with i = 1. There is a stable fixed point at 0.5173 
with /'(0.5173) = 0.2973 so that (Dl) is satisfied with k = 1, and clearly 
(D2) is also satisfied. Therefore, by Theorem 3 the limiting eigenvalues are 
{0.2973" : n > 0} = {1, 0.2973, 0.0884, 0.0263, . . .}. 

The next four examples correspond closely to the values considered by 
Tateno and Jimbo. 

Example 3. / = 2 and k = 0.2. Figure 5 shows / and There is a 
period 2 stable orbit {0.3527,0.7593} with (/2)'(0.3527) = 0.7445 and a 
period 2 unstable orbit {0.4654,0.9329} with (/2)'(0.4654) = 1.5043. By 
Remark 3 following Theorem 2 the limiting eigenvalues are {(0.7445")^/^ : 
n > 0} U {(1.5043""~i)^/2 . ^ > q} ^ ±0.8628, ±0.8153, ±0.7445, . . .}. 

Example 4. / = 2 and k = 0.5. Again we show / and see Figure 6. There 
is a point of discontinuity with D = {0.5489} with E = {0.8567,0.3057}. In 
this example f~^{D) = {.1174} / 0, but f-^{D) = /"^({0.1174}) = 0, and 
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(D3) holds with 1 = 2. There is a period 2 stable orbit {0.3651,0.6586}, 
and the product of the derivatives along the orbit is (/^)'(0.3651) = 0.2544 
so that (Dl) holds with k = 2. It is clear from the plot of that (D2) also 
holds. Thus Theorem 2 implies the limiting eigenvalues are {±(0.2554")^/^ : 
n > 0} = {±1, ±0.5044, ±0.2544, . . .}. This is one case considered in Figure 
4 of Tateno and Jimbo [16] and our predicted limiting values can be seen at 
the extreme left edge of the k = 0.5 part of Figure 4 of [16]. 

Example 5. / = 2 and k = 0.8. Figure 7 shows / and f^. We are again in 
the setting of Theorem 3. There is an attracting period 3 orbit 

{0.4218,0.6330,0.7352} 

and the product of the derivatives along the orbit is (/2)'(0.4218) = 0.088076 
so the hmiting eigenvalues are { (0.088076" . > q} ^ {u;''0.4449" : r = 
0, 1,2 and n > 0} where uj = e^'^'^l'^ is a cube root of unity. These can be 
seen at the extreme left edge of the k = 0.8 part of Figure 4 of [16], and also 
in Figure 5(c) of [16]. 

Example 6. / = 2 and k = 0.9. Figure 8 shows / and There is now an 
attracting period four orbit {0.4378, 0.6236, 0.6978, 0.7480} and the product 
of the derivatives along the orbit is 0.043991 and the conditions of Theorem 
3 are satisfied so the limiting eigenvalues are { (0.043991'' : n > 0} = 
{^'"0.4580"^ : r = 0,1,2,3 and n > 0}. This gives a theoretical justification 
for the k = 0.9 part of Figure 5(c) of [16] where fourth roots of one appear 
in the limiting spectrum. 

6. Proofs for Section 2. We show first that it suffices to prove the 
results for the case of constant input function. For any constant / define 

k{t) = e-^^ l\^'{I{s)- I)ds. 
Jo 

Then Xf = Xf - k{t) satisfies 

dX^ = (^-jXf + I^dt + edWt. 

Define g{t) = g[t) — k{t), then started at xq — k{tQ) at time to hits the 
threshold g{t) at the same moment that Xf started at xq at time to hits 
the threshold g{t). Moreover if ^, /, m, Q, a, and P are defined using the 
constant input I[t) = / and the threshold function g in the same way as ^, 
/, m, Q, a, and P are defined using the original input I{t) and original 
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threshold then ^(t|to, a^o - ^(^o)) = §it\to, xq) - k{t), f{tQ,XQ-k{tQ)) = 
/(to, xq), m{to, xo - k{to)) = m{to,xo), § = {{to.XQ - k{to)) : {to,xo) G Q}, 
a{to,xo - k{to)) = a{to,xo) and pto,^o-fc{to)(fe (z A) = p*o,^o(^^ ^ a) for 
any Borel subset A C K. Therefore any of the results of Section 2 proved 
under the assumption that I{t) = I can be converted into the corresponding 
result for a more general function I{t). 

This method of converting a problem with a time varying input into 
one with a constant input and time varying threshold (and reset) is used 
in Scharstein [12]. For the reminder of this section we shall assume that 
I(t) = I, so that is either an Ornstein-Uhlenbeck process (if 7 > or 
else a Brownian motion with constant drift. 

Proof of Proposition 1. Suppose that {to,xo) G Q. An application of 
the implicit function theorem to F{s,x,t) := (,{t\s,x) — g{t) at the point 
(to, a^O) C(/(^0) a^o))) implies that f{s,x) < 00 in some neighborhood U of 
(to,XQ) and that £ C^(f/, M). The continuity of g implies that {{s,x) : 
X < g{t)} is a neighborhood of (to, 2:0), ^-^d the fact that g is and / is 
continuous on U implies that {{s,x) : g'{f{s,x)) — I + jg{f{s, x)) < 0} is a 
neighborhood of (to, Xo). □ 

Proof of Proposition 2. The compactness of G implies the existence of 
61 > and 82 G (0, 6] such that 

C(t|to, Xq) < g{t) - 5i for to < t < /(to, xo) - 5 

and 

^(t|to, Xo) > g{t) + 61 for t = /(to, Xo) + 52 
for all (to,xo) G G. These two inequalities imply that 

P*»'"°(|T"-/(i0,X0)| ><^) 

< P*0'"° {{XI > g{s) for some s G [to, /(to,xo) - 6)} 
U {X(f(to,xo)+s,) < 9ifito,xo) + 62)}) 

<F^«'*o( sup |Xf-e(s|to,xo)| ><5i 

\to<S<f{to,Xo)+S2 

From (4) and the definition of C(s|to,xo) we have 

\X! - C(s|to,xo)| < 7 r \X'{u) - au\to,xo)\du + e\Ws - Wt,\ 
J to 
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and then Gronwall's inequality gives 

sup |Xf-e(s|to,^o)| <£e^^*"*"^ sup \Ws-Wt,\. 

to<S<t to<S<t 

Define T = sup{f{tQ,XQ) + 52 — to : {to,^o) £ G}, and note that the com- 
pactness of G implies that T < oo. We have 

sup \XI - ^{s\to,xo)\ < ee''^ sup 

to<S<f{to,Xo)+S2 to<S<to+T 



dist 



ee^'^ sup I Ws 

0<s<T 



Therefore 



\to<S<f(to,Xo)+52 J 

< P (ee^^ sup \W,\ > 6i] 

\ 0<s<T J 

< 2P (^Wt\ > ^^^^ 

and the result follows directly. □ 

The proof of Theorem 1 is rather lengthy so we will split the remainder 
of this section up into several subsections highlighting the main components 
which will be tied together in Section 6.5. We begin with a simple example 
for motivational purposes. 

6.1. Non-leaky case with constant threshold. Suppose that 7 = and 
/ > and that the threshold g{t) = B is constant. Then is just Brownian 
motion with drift / and G = {{s,x) : x < B} with f{t) = to + (B — x)/I. 
It can be shown (see for instance [1]) that for any xq < B the first-passage 
time density is given explicitly by 

p%t\to,Xo) = --^q'{t\to,Xo) 

t — to 

for t > to where q^{t\to,xo) is the transition density at B for the random 
variable Xf given X^^ = xq. Thus 

euu \ ^-a;o [ {B-xo-I{t-to)y\ 

p%t\to,xo) = ^ , — ^exp' ^ 



^£(t_to)3/2 ^1 2E\t-to) 

I{f{to,xo)-to)) i P{t-f{to,xo)) 



V2^e{t - to)3/2 ■ [ 2e2(t -to) j ' 
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Replacing t by f{to,xo) + et we get 
ep%et + f{to,xo)\to,xo) =- 



27r{f{to,xo)-to + et)^/^ 
X exp < — 



2{f{to,xo)-to + et) 



I{f{to,xo)-to)) j IH^ 

exp ' 



2^(/(to, xo) - to)3/2 " i 2(/(to, xo) - to) j 
= p^(t|to,2;o) 

as e — > because in this setting a'^{to, xq) = f{to, xq) — to and m{to, xo) = I. 
This shows the pointwise convergence imphed by (6). The full strength of 
(6) will follow from the techniques developed below for the general case. 

6.2. Durbin's Theorem. If 7 > or g is not constant, then we no longer 
have explicit formulas at our disposal. However, since is a Gaussian 
process, we have the following result of Durbin [4, page 100], valid for any 
7 > and any function g. 

Theorem 4. Suppose that xq < g{to). For t > to 
(13) p%t\to,xo) =b%t\to,xo)q'{t\to,xo) 

where q^{t\to,XQ) is the density under of Xf evaluated at g{t), and 

b'it\to,xo) := lim-^E*<"-°[l,.>,(<7(s) = 5(t)]. 

s/t t — S 

We call the density term and b^ the slope term in the decomposition 
(13). 

6.3. Analysis of the density term. We deal with q^ in much the same 
way as we dealt with p^ in the example from Section 6.1. Define 

Jto [ t - to if 7 = 0. 

Then 

^''^ ' = V2^eait\to) ^"^^ 2.V^(t|to) ^ 

where (,{t\to,xo) is the solution to the ODE x' = — 7X + I for t > to with 
S,{to\to,xo) = Xo- Note that a'^(to,xo) = a'^{f{to,xo)\to)/m'^{to,xo). 
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Lemma 1. If G is a compact subset ofQ then there exist 6,eo,K,ai > 
so that 

(15) \m{to,xo)£q%f{to,xo) + u\to,xo) -pr{u/e\to,xo)\ < Kee'^''/^''''^? 
and 

(16) (e + \u\)q'{f{to, xo) + n|to, xq) < i^e-^'/^eV? 
for all (to, xq) £ G, \u\ < 5 and e < Eq. 

Proof. The compactness of G implies the existence of (5i > such that 
f{to,xo) > to + 2di for all {tQ,xo) E G. Noting that 

9{f{to,xo)) - C(/(to, xo)\to,xo) = 

and 

g'ifito, Xo)) - C'ifito, xo)\to, Xo) = -m{to, xo) 
we obtain by Taylor's theorem 

[9{f(to,xo) + u)-£,{f{to,xo) + u\to,xo)f , , 

27777 — u~\ = ~Ul 7 + uKi[to,xo,u)) 

a'^[f[to,xo)+u\to) a^{to,xo) 

where the remainder term satisfies 

|i?i(to, 3^0, ^^)| < Ki for all {to,xo) G G and |n| < 5i 

for some Ki. Similarly we have 

m{to,xo) 1 /I , D ^^ 

-TTTT r- iTT = — r. V (1 + -"^2(^0, xo,u)) 

cr{f{to,xo)+u\to) ar{to,xo) 

where the remainder term satisfies 

\R2(to,xo,u)\ < K2 for all {to.,xo) G G and |n| < 5i 
for some K2. For ease of notation in the following calculation we drop the 
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arguments of cTj- and Ri and i?2- We have 
\m{to, xo)eq'^{f{to, xq) + u\to, xq) - pr{u/e\to,xo)\ 

(1 + tx/?2)e-"'('+"''^)/'"'"' - e-"'/2.2.? 



27r(7T 



< 



27r(7T 



+ 



max {e-^(i+««i)/2-=-? , e-V2.^^? } j 



At this point choose 6 < 6i so that -fCif^ < 1/3. Then for |u| < 6 we have 
|m(to, 2;o)eg^(/(to, a^o) + '"1*0, xq) - Pr{u/e\tQ, xq)\ 



< 



1 



27r 



K, 



u 




u 




ear 


2 


ear 





< 



K2 + ^]K,e-y'^''^'^ 



where 



K3 = max ^max{|x|e ^ ^^"^ : x £ M},max{|xpe 



The resuh (15) follows directly, with af = 2 max{cj^(to, a;o) : {to,xo) G G}. 
The proof of (16) is similar (and simpler) and is left to the reader. □ 

6.4. Analysis of the slope term. The slope term 

b%t\to,xo) := hm ^E*"'-o[i^,^^(^(,) _ xi)\X^ = g{t)] 

s/t t — S 

involves the pinned process given = xq and = y. For t^ < t define 

sinh7(s — to) 



(17) 



ip{s\to,t) 



sinh7(t — to) 

s-tp 

t-to 



if 7 > 
if 7 = 
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and 

(18) n{s\to,xo, t, y) = C{s\to, xo) + ^(s|to, t){y - ^{t\to,xo)) 

for tQ < s < t. Our next proposition gives us a useful representation for the 
pinned process. 

Proposition 4. For to < t the conditional distribution of {X^ : to < 
s < t} given = xq and Xf = y can be written 

X^ = fi{s\to,xo,t,y) +£Us, s£[to,t] 

where K{Us) = and 



E( sup m]<i ^\/(^^^^""^ - 1)/^^ 

Kio<s<t J I K{t-to) if 7 = 0. 



The proof of Proposition 4 rehes on the following standard result regard- 
ing the conditioned law of a Gaussian process. We include the proof for 
completeness. 

Lemma 2. Suppose that {Zg : s £ S} is a real valued Gaussian process 
defined on some set S CU. with mean and covariance p{s, t). Let N > 1 
and ti,t2, ■ ■ - tN S -S*, and suppose that the matrix A := {p{ti,tj))ij=i^,,,^j\f 
is invertible with = B. Then the conditional law of {Zs : s G 5} given 
Zt^ = zi, Ztj^ = ZN is given by 

N 

+ Pi^^ti)Bij{zj - H{tj)) + Us, s £ S, 

where 

N 

Us = {Zs - /i(s)) - p{s,ti)Bij{Ztj - fi{tj)). 

The process {Ut ■ s £ S} is Gaussian with mean and covariance function 
p{s,t) = p{s,t) - Y.f,j=i p{s,ti)Bijp{t,tj). 

Proof. Define the process 

N 

Vs = p{s) + Y p{s,ti)Bij{Ztj - p{tj)) 
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for s £ S. Clearly {Vs : s G 5} is measurable with respect to a{Zt^, ..,Ztj^} 
and is a Gaussian process. Moreover Vt- = Zt- for j = 1,2,...,N. Now 
define Us = Zg — Vg. The process {Us : s E 5} is a mean zero Gaussian 
process, and a direct algebraic calculation gives that Cov{Us, Zt-) = for 
all j = 1,2, . . . , N , so that {Us : s G 5} is independent of a{Zt^, .., Ztj^}. 
Therefore the law of Zs = Us + Vs given Zt- = Zi, for z = 1, is the same 
as the law of Us + /^(s) + J2f,j=i Pi^^i U)Bij{zj — fi{tt), and the first assertion 
is proved. The covariance of {Us : s £ S} is a direct algebraic calculation 
and is left to the reader. □ 

Proof of Proposition 4. We will apply Lemma 2 to the process {X^ : s > 
to} given by 

dXI = {-jXI + I)ds + edWs 

with initial condition Xf^^ = xq, conditioned by Xf = y. We will take = 1 
and condition at the single point t, so that = 1 and Bu = p(t,t)~^. For 
s > to we have fi{s) = ^{s\tQ, xq), and for si,S2 > to we have 



P{S1,S2) 



(eV27) (e-^l^'i-'^l - e-^^'^i+^^-^to) j if ^ > o 
e"^ min(si — to) •S2 — ^o) if 7 = 0. 



Lemma 2 implies that the conditional law of {X^ ■ to < s < t} given X^ 



IS 



Ks\to,xo,t,y) + Ul 



where 



and 



n{s\to,xo,t,y) = £,is\to,xo) + 



p{s,t) 



Ul = Xl-Cis\to,xo) 



p{t,t 
p{s,t) 



{y - Cit\xo,to) 



p{t,t) 



iXl-mto.xo)). 



Case 1: 7 = 0. We have X^ — ^{s\tQ,XQ) = e{Ws — Wto) and p{si,S2) 
min(si — tQ,S2 — t^), so that C/| = eUs where 



Us 



W,-Wt 



to 



to 



t-to 



{Wt-Wt,) 



The process {Us : to < s < t} is the Brownian bridge on the interval [to)i]- 
For < u < 1 define Wu = (Wt^+uit-to) ~ ^to)/V* ~ then W is standard 
Brownian motion and 



Us = 



Wt 



(s-to)/{t-to) 



to 



Wi 



-t-to, 

imsart-aap ver. 2009/05/21 file: IF_Article.tex date: July 22, 2009 



20 BAXENDALE AND MAYBERRY 

We have 

sup \U{s)\ = ^yt — to sup IVF^ — uWil 

to<s<t 0<u<l 

and in particular 

El sup \U{s)\] =Vt^Ei sup |W^n-itM^i|) =K^/t^. 

\to<s<t I \0<'U<1 I 



Case 2: 7 > 0. For < s <t -we have 

Jto 

and 

p{t,t) ~ (1 _ e-27(t-to)) 
Therefore Ug = eUg where 



and so 



'~Jto ^ (1 _ e-27(t-to)) X/ 



For fixed tQ < t define a : [tQ,t] ^ M hy 

""^^^ " (1 _ e-27{t-to)) • 

The function a is continuous and strictly increasing from [t(),t] onto [0,1], 
with inverse 

a-\u) =to + ^ log ((1 -u)+ ne2^(*-*o)) , < u < 1. 
For < li < 1 define 
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Then {VK(n) : < n < 1} is a standard Brownian motion, and 



I g27i _ g27to 



27 



{w{a{s))-a{s)W{l)). 



Therefore 



sup \U{s)\ < e-^*o sup |eT"C/(s) 

to<s<t to<s<t 



/g27(t-to) _ I 



sup 

27 to<s<t 



sup 

^7 o<M<i 



Ty(a(s)) -a(s)W(l) 
W{u)-uW{l) 



In particular we have 



E sup \U{s)\ < Ke^ 

\to<S<t J 



o2j{t-to) _ I 



27 



□ 



For to < t define 

(19) f3{s\to,xo,t) = g{s) - fi{s\to,xo,t,g{t)) to<s<t, 

so that (3{s\tQ,XQ,t,) is the gap between the threshold g{s) and the con- 
ditional expected value of given Xf^ = xq and X^ = g{t). A direct 
calculation gives 

P'{s\to,im) = g'{s)+-fC{s\to,xo)-I-ip'{s\to,t){g{t)-^{t\to,xo) 



9(t)-xo '^ ^ ^ g{t) - xo \^ ^ _^ -f{t -to) cosh (t - s) 



(21) 

if 7 > and 
(22) 

if 7 = 0. 



+ {ig{t)-l) 



t — to J V t — to J \ sinh7(t — to) 

/ cosh7(t — s) — cosh7(s — to)" 



l3'{s\to,xo,t) = {g'{s) 



sinh7(t — to) 

_ gjt) - xq 
t-to 
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Following Durbin [4] we can write 
(23) b%t\to,xo) = bl{t\to,xo) - b%t\to,xo) 

where 



and 



blit\to,xo) := lim-^E*0'-«[(5(s) - = git)] 

s/t t — S 



¥{t\to,xo) :=lim(t-5)-iE*«'^«[l,.<,(5(s)-Xf)|Xf =9(t)]. 

s /"t 



By Proposition 4 we have 

Oi(tto,a;o) = hm — 

s/t t — S 

= -(3'{t\to,xo,t) 

(24) = -g'{t)--f^{t\to,xo) + I-^p'{t\to,t){g{t)-^{t\to,xo)). 

In particular bi{t\to,xo) is independent of e, and henceforth we shall write 
bi{t\tQ,xo) = bi{t\to,xo). Notice that 

(25) 6i(/(to, xo)\to, xo) = -g'ifito, xo)) - 75(/(to, xq)) + I = m{tQ, xq). 

The following result of Durbin [4] will enable us to control the W term in 
(23). 

Proposition 5. For any t > to we have 

(26) b'{t\to,xo) = f\i{t\r,g{r))p'{r\to,xo,t,g{t))dr. 

where p'^{r\to,xo,t,y) denotes the conditional density function for t'^ given 
that Xf = y and X^^ = xq 

Proof: This a corrected restatement of [4, Section 6, equation (27)]. The 
Markov property for X^ gives 

b%t\to,xo) = lim Te-'^W f 9is) -^"(') X^it)=g{t))f{r\to,xo,t,g{t))dr 
s/t Jto \ t- s J 

^ r fmri9MA)pe(^r\to,xo,t,git))dr. 

s/t Jto \ t - s J 

It is easily checked from the expressions (21) and (22) that \l3'{s\r, g{r),t)\ 
is bounded for to < r < s < t, and since we also have (3{t\r, g{r),t) = 0, the 
passage of limsyt inside the integral is justified by the bounded convergence 
theorem. □ 
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Remark 4. The substitution of (26) into (23) gives an integral expres- 
sion for b^{t\tQ, xq) in terms of the conditional density function j3^(r| to, xo,t,y) 
The relation (6.1) then gives an integral equation for the first passage den- 
sity p'^{t\tQ,XQ) which is a special case of the integral equation of Buonocore, 
Nobile and Ricciardi [1]. In particular [1] can be used to give an alternative 
derivation of (24). The paper [1] deals with time-homogeneous diffusion pro- 
cesses, and the paper [4] deals with Gaussian processes, and we are working 
in the intersection of these two classes of processes. 

Proposition 6. Let G be a compact subset of Q. Then there exist 5 > 
and K < oo such that 



whenever (to, xq) G G and \t — /(to, xo)\ < 6. 

We prepare for the proof with the following pair of Lemmas. 

Lemma 3. Given the compact set G C G there are positive 6, 6i, ki and 
k2 such that 



whenever (to, 2;o) £ G and |t — /(to, xq)\ < 5. 

Proof. We give the proof for the case 7 > 0; the proof for the case 7 = is 
essentially the same. Equation (21) shows that (3'{s\to,xo,t) is a continuous 
function of {tQ,XQ,t, s) on the set where t > to. Also, putting s = t = 
/(to,xo) in (20) gives 

/5'(/(to, xo)\to,xo, /(to, xq)) = g'ifito, xo))+7fl'(/(io, xo))-I = -m(to, xq) < 

for (to, 2:0) G G- The compactness of G gives 82 > and ki > such that 
/3'(s|to, xo, t) < —ki whenever (to,xo) G G and f{to,XQ) — 62 < s < t < 
f{tQ,xo) + 52- Since /3(t|to, xo, t) = it follows that 



whenever (to,xo) G G and /(to,a;o) — < s < t < f{to,xo) + 62- 

The definition of f{to,XQ) as the time of first intersection of ^(t|to,xo) 
with g{t) implies that 



(27) 



(t|to,xo) - 6i(t|to,xo)| < eK 




ift-5i<s<t 
if to < s < t — 61 



(28) 



P{s\to,xo,t) > ki{t - s) 



min{/?(s|to, xo, /(to, xq)) : to < s < /(to, xq) - 62} > 
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for all (tQ^xo) G G. The compactness of G gives the existence of iJa > and 
k2 > such that 



(29) 



P{s\to,xo,t) > k2 



whenever {to,xo) E G and \t — f{tQ,xo)\ < 83 and to < s < f{to,xo) — 62- The 
result is now a simple consequence of (28) and (29), with 6 = min(J2/2, (^3), 
5i = 62/2. □ 

Lemma 4. Let H be a compact subset o/M. There is Ki such that 

\bi{t\r,g{r))\ < Ki{t - r) 
whenever r,t £ H and r < t. 

Proof. Again we give the proof for the case 7 > and leave 7 = to the 
reader. Putting s = t and {tQ,xo) = {r,g{r)) in (21) gives 



bi{t\r,g{r)) 



9it) - 9{r) 



9'{t) + 



+ 7 



t - r 

/ cosh[7(t — r)] — 1 
V sinh[7(t — r)] 



7(^ - r) 
sinh[7(t — r)] 

{9{t)-I)- 



1 



g(0 - 9{r) 
t - r 



Using the inequalities 
cosh u — 1 



sinhn 
we get 

\hi{t\r,g{r))\ 
t — r 



< — and 
6 



u 








< min 




sinhn 







< 



g{t) - g{r) - [t - r)g'{t) 



+ 



7 

V6 



\u\ 



9{t) - 9{r) 



7 



+ylffW-/| 

1/7 / 7^ 

< 7T sup |/(s)| + sup \g'{s)\ + —\g{t) - I\ 

^ r<s<t VD r<s<t ^ 



and the result follows from the compactness of H. 



□ 



Proof of Proposition 6. Let (5, 5i, ki and k2 be as in Lemma 3, and then 
Ki as in Lemma 4 with H = {t : to < t < f{tQ, X())+5 for some (to, xq) G G}. 
Now fix (to,a^o) £ Gandt G [f{tQ,xo) — 5,f{tQ,XQ)+5] and apply Proposition 
4 with y = g{t). Under the conditions X^^ = xq and Xf = g{t) we have 

= inf{s > to : = [i{s\tQ,XQ,t)}. 
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Define \\U\\ := sup(g<^<JJ7s|- If ^11^^11 = x < min(/ci5i, ^2) then by Lemma 

3 we have eUs < (3{s\tQ, xo,t) for all s < t — x/ki and so > t — x/ki and 
|5l(^|T^c^(T^))| < Kix/h = eKi\\U\\/ki. li e\\U\\ > m.m{ki6i,k2) we can 
use the estimate \bi{t\T^ , g{T^))\ < Ki{t — t^) < K{t — to) from Lemma 4. 
Together we have 

61 t r^5 < e U max —, • 

\ki mm{kidi,k2) J 

Proposition 5 gives 

rt 

\b'{t\to,xo)\ < / \bi{t\r,g{r))\p'{r\to,xo,t,g{t))dt 
Jto 

= E*0'-o(|6l(^|r^ff(r^))||Xf =5(t)) 
Ki Kljt -to) \ 
ki ' min(/ci(5i, ^2) / ' 

and the result now follows by Proposition 4 and the compactness of G. □ 

6.5. Proof of Theorem 1. We are now ready to complete the proof of 
Theorem 1. Given G, choose (5 > sufficiently small and K < 00 sufficiently 
large so that the results of Lemma 1 and Proposition 6 are valid. By Theorem 

4 we have 

k/(/(io,2;o) +u\to,xo) - Pr{u/e\to,xo)\ 
< \m{to,xo)eq%f{to,xo) + u) - pr{u/e\to,xo)\ 

+ \b%f{to,xo) + u\to,xo) - m{to,xo)\eq^{f{to,xo) + u\to,xo) 
= I + 11. 

Now by Proposition 6 

\b%f{to,xo) +u\to,xo) - m{to,xo)\ 

< \b^{f{to,xo) + u\to,xo) - bi{f{to,xo) + u\to,xo)\ 
+ \h{f{to,xo) +u\to,xo) - bi{f{to,xo)\to,xo)\ 

< eK + \u\Ki 

where 

Ki = sup{\b[{s\to,xo)\ : |s - f{to,xo)\ < 5} < 00. 

(The finiteness of Ki uses the fact that f{to, xq) — to is bounded away from 
on the compact set G.) The result now follows from Lemma 1, using (15) 
on / and (16) on //. □ 
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7. Proofs for Section 3. In what follows, d denotes the standard quo- 
tient metric on S induced by the Euclidean metric on M. The starting point 
for the proof of Theorem 2 is the following splitting of S. 

Proposition 7. There exist neighborhoods Vi := Bs^{6u), V3 := Bs^Og), 
and constants 5 > 0, iV G N such that 

(1) d{l{e), Vi) > 6 for every 9 Vi 

(2) dif le), Kf) > 5 for every 6 £ V3 

(3) For every 9 e V2 := E>/{Vi U V3), we have P{9) G V3, Vn > N. 



Proof. This is the same as Proposition 1 in [10]. □ 

We can write any (p E -B(S) in the form 4> = (t>i+ (t>2 + (f>3 where = (j)lvi ■ 
The action of can then be described by the block decomposition 



^11 


rpe 


rpe 

^13 


rpe 


rpe 
^22 


rpe 

^23 


rpe 


rpe 


rpe 

^33 



where 

(30) Tt^m) = J mpuomde 

with plj{9\eo) = lv,{eo)lv^{9)p^{9\eo). The choice of the sets Vi, V2 and ^3 
implies that 



T?i 


T?2 







rpO 

^22 


rpO 

^23 








rpe 
-^33 



For e > write 














rpe 


rpe 

^12 


rpe 

^13 


rpe 
^21 








and Tip = 





rpe 
^22 


rpe 
-'23 


rpe 

^31 


rpe 

^32 













rpe 

^33 



Lemma 5. There are finite positive constants K andM such that ||Tj^ ||oo < 
Mee-^/"^ for ij = 21 or 31 or 32, and ||(T|2)^+i||oo < Mee"^/^". 

Proof. For ij = 21 or 31 or 32 we have by Proposition 7 

||7;^J = supP'^(Gf G Vj) < supF^{d{eiJi9)) > 6) 

< sup F\\T'-fit)\>5), 

0<i<l 



imsart-aap ver. 2009/05/21 file: IF_Article.tex date: July 22, 2009 



SPECTRAL ANALYSIS OF SIFS 



27 



and the first set of results follows by Proposition 2. Also by Proposition 7 

ll(r|2)^+'||oo < sup P^(d(e^^+i,/™(0)) > 6), 

and the second result follows using Proposition 2 together with the inequality 

N 

i=o 

where L = sup|/'(6l)|. □ 

It follows directly from Lemma 5 that for any r > there is eo > 
such that llr^plloo < r and all eigenvalues of have modulus less than r 
whenever e < eo . In order to complete the proof of Theorem 2 it suffices to 
describe the eigenvalues of the operators Tf^ and Tgg as e — 0. This will be 
carried out in Sections 7.1 and 7.2. 

7.1. Behavior near a stable fixed point. Recall that is the restriction 
of the transition operator to a neighborhood V3 of the stable fixed point 
6s, and that f'{6s) = Cg- The main result of this section is 

Proposition 8. Every non-zero eigenvalue of T^^ is of the form c" + 
0(e) as e ^ for some n >0. 

We can reparameterize § so that = and V3 = {—6s,Ss)- Then /(O) = 
and /(O) = Us for some Us > 0. Also /'(O) = /'(O) = c^. Recall that Tgg is 
the operator on BiVs) defined by 



(31) Tls(l){to)= J lv,it)fii\to)Ht)dt, toGVa 
where 

(32) f{t\to) = Y.f{t + n\to) 

We then extend T33 to an operator on B(M.) via 

(33) Tl^^ito) = lv,ito) f lv,it)p'it\to)mdt, toe 



and look at the re-scaled version := (Ue) ^ o T33 o Ue where Us4>{x) = 
(j){x/e), so that 

rj0(to) = Iviito) J lvg{t)HiW{et\eto)dt 
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with = Vi/£ = {—5s/s, Ss/e)- For ease of notation, we drop the subscript 
s from Cs and 5s and the subscript 3 from V3 and V^. 
From (32) we get 



£p'^{et\£tQ) = ^ ep^ {et + n\£to) . 

Let 

Pmmn(i|*o) = e/(ei + n^leto) 

denote the main term in this sum. Theorem 1, together with the hmiting 
behavior m(eto) — ^ m(0) and cr(eto) — > o'(O) and (/(eto) — ''^s)/^ — > cIq 
imphes that 

(34) fma^nm) - Vr{t " Cto|0) = ^ -{t-ct,)V2.li^) 

V27r(j^(0) 

in some sense. (For details of this calculation see equation (35) later.) This 
suggests that in some sense, — > Tg as e — > where Tg is the operator with 
kernel pT-{t — cto\Q). 

In order to make this precise, for /c G M define \\(t)\\k = sup{|(/>(x)|e~'^^ : 
X € M} and Wk = {(^ : M — > M : is measurable and ||</>||fc < 00}. Then 
Wk with the norm || • \\k is a Banach space. Let CiyVk) denote the set of all 
bounded linear operators T on Wk with operator norm 

llrllfe = sup{||T</.||fc -.c^eWk and < !}• 

Proposition 9. For all k > sufficiently small, we have = Ts+0{£) 
inC{Wk). 

Proof. Define the operators T^aj^, T^^t with the following kernels: 
T^^, ^ lv4t)lv4to)pr{t - Cto\0) 

and write 

rp£ _rp _ frpe _ rpS \ . (rpE _ T^^ ^ _|_ (J^l^ — T ^ 

s s y-^ s ^mamJ ' K-^mam -^cut) ' \^ cut ^s) 

= I + 11 + III 

The proof will consist of bounding /, //, /// as operators on Wk for small 
enough k. 
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To bound /, we note that for (p G Wk we have 



\iT!-T^a^n)m)\ 



J lvieto)lvi£t)(l){t)ep''{et + n\eto)dt 
< Uhe''"^'hv{eto) I J2 J lv{et)ep%et + n\eto)dt 
= Uhe''"^'hv{eto)P'''^ It' e [j V + n 



Now |eto| < 5 imphes |/(eto) — ns\ < 5, so that 



e U V + nW {\t' - /(rfo)l > 1 - 25}. 

n^ns ) 

We can assume without loss of generahty that 5 was chosen small enough 
so that 5i := 1 — 25 > 0. Then Proposition 2 gives 

\{T! - T^aMt^)\ < Uhe^'^'^'P^'^ (|r^ - /(eto)| > 1 - 25) 
and so 

This is at most 0(e) as long as /c < K^^. 

The calculation giving an 0{e) bound for /// concerns the effect of the 
cutoffs lye(to)lyE(i) on a Gaussian kernel. This is the same as in the proof 
of Equation (17) in [10] and is omitted. 

For // we have 

P"ma^n{t\to) = ^/(/(O) + rf|rfo) 

= ep'Uieto) + et- {f{eto) - /(0))|£to) 
= ep%f{eto) + et- f{eto)\eto) 

= Pr{t-cto\0) 

+ [ef{e{t - f{to)) + f{eto)\eto)-Pr{t - r{to)\eto)] 
+ [pAt - r{to)\eto) - Pr{t - ctolO)] 
(35) = pr{t - cto\0) + rl{t\to) + rUt\to) 
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where /^(io) = e~^/(£io)- Define operators T^. with kernels 

ly.(t)ly.(to)rf(t|to), 

i = 1,2. The operator T^^ deals with the effect of changing the mean /^(to) 
and standard deviation ar{£to) of a Gaussian kernel, and Equation (14) in 
[10] gives an 0{e) bound for HT^jUfc for all sufficiently small A; > 0. For T^_^, 
suppose A; > and (j) £ Wk- We have 

\Tr'Mto)\ = J lv4t)^v4toVi{t\to)<t>{t) dt 
J ly.(t)ly.(to) 

X [ep%e{t - f{to)) + f{eto)\eto) - Pr{t - f {to)\eto)](j){t) dt 

In order to apply Theorem 1 we need to restrict to \et — f{eto)\ sufficiently 
small. This can be achieved for tQ,t G by choosing 5s sufficiently small 
in Proposition 7. Then 



K£\(\)%\yeiytQ) =exp^ 



1 - 2ka\ I 1 - 2^^^! J 

Furthermore by shrinking V if necessary, we can find ci < 1 such that 
I /''(to) I < ci|to| for to e ^^ and then 

< KeUUly.ito) ^^"^ ^ exp ^ ^ 



It follows that 



11 V 27r(Ti 



1 - 2/ccjf 



so long as cf < 1 — Ikaf, that is, A; < (1 — cf )/2crf . This completes the proof 
of Proposition 9. □ 

Now Ts is the transition operator for the Markov chain X„ = cX„,_i + Xn, 
and has eigenvalues c", n > 0. It follows from Proposition 9 together with 
standard perturbation results for linear operators (see Kato [7]) that the 
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eigenvalues of acting on Wk are of the form c" + 0(e). Finally, since the 
operator involves the indicator functions lye (to)!^^ (i) the eigenvalues 
and eigenfunctions do not depend on the value of k, and in particular the 
non-zero eigenvalues of acting on Wk coincide with the non-zero eigen- 
values of T^^ acting on B{V). This completes the proof of Proposition 8. 

7.2. Behavior near an unstable fixed point. Here Tf is the restriction of 
the transition operator to a neighborhood Vi of the unstable fixed point 
6u, and f'{Ou) = Cu with \cu\ > 1- The main result of this section is 

Proposition 10. Every non-zero eigenvalue ofTf-^ is of the form 



as e — > for some n > 0. 

Here we localize in the neighborhood Vi of the unstable fixed point and 
replace and Tg in Section 7.1 with and defined in the same way, but 
using Vi in place of V3. We obtain a perturbation result similar to Propo- 
sition 9, but in a class of spaces of functions with exponential decay. Note 
that all calculations in the previous section before the proof of Proposition 
9 are valid if |c| > 1 as well. 

Proposition 11. For all k > sufficiently small, we have = Tu + 
0{e) inC{W.k)- 

Proof. Define the operators T^^j^ and T^^^ as in the proof of Proposition 
9 and decompose — in the same way as / + // + ///. The bound on 
/ is obtained in essentially the same as in Section 7.1; the only difference 
is that now etQ G Vi implies |/(eto) — n^l < 5 for some 5 > 0, so that we 
need 5i := 1 — 5 — 5 > 0. The bounds on /// and the bound on T^^ in 
the decomposition of // involve Gaussian kernels and are the same as those 
used in the proof of Theorem 9 in [10]. Finally, by shrinking Vi if necessary, 
we can find ci > 1 such that \ f^{to)\ > cito for all eto G Vi. By a similar 
argument to that given in Section 7.1 above, for k > and (j) E W-k we 
obtain 



-1 



+ 0(6) 



T^\(l>{to)\<KeU\Ulv^{to) 



+ 2kal 



exp 



( 



1 + 2A;o-2 




) 



It follows that 
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SO long as cf > 1 + 2kaf, that is, A; < (cf - l)/2a^. □ 

Since has eigenvalues |c|~-'^c~" if |c| > 1 (see Section 5 of [10]), the proof 
of Proposition 10 again follows from standard perturbation arguments. 

8. Proofs for Section 4. Proof of Proposition 3. Since D is finite 
it suffices to prove the result separately for each Iq G D. Suppose first 
that Iq & D satisfies (C')(i), so that /*(to) = /(^o)- Since ^ is a continuous 
function of both arguments and g is continuous at s = tj, there exist 6,6i > 
and 82 € (0, 6] such that 

as\t) < g{s) -61 for t<s< f{to) - S 

and 

^{s\t) > g{s) + 61 for s = /(to) + 52. 

whenever |t — to| < ^- As in the proof of Proposition 2, these two inequalities 
imply that 

r\\T'-f{to)\>5) < pM sup |X^(s)-e(s|0| ><5il . 

\t<S<f{to)+52 J 

If instead to £ D satisfies (C')(ii), then again using the continuity of ^ and 
g, there exist 6 > 0, 61 > and 62 G (0, 6] such that 

C{s\t) < g{s) - 61 for s G [t, /(to) - <5) U (/(to) + 6, /*(to) - S) 

and 

> 9{s) + 5i for s = /*(to) + <52. 
whenever |t — to| < ^- These two inequalities imply that 

IP*(rf(T^{/(^o),r(^o)})><5) < P*(K</(to)-5} 

U{f{to) + 6<T' <r{to)-6} 

uK>/*(to) + M) 

< I'M sup |X^(s) -e(s|t)| > (5i ) . 

\t<s</*(to)+<52 / 

In either case the proof is completed using the same method as in the proof 
of Proposition 2. □ 

The following Lemma will be used several times in the proof of Theorem 
3. Recall that / = / mod 1. 
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Lemma 6. Assume that f satisfies (A), (B') and (C). For 9 £ S 
suppose f^{6) D for < i < n — 1 and let V be any neighborhood o//'"(0). 
There is a neighborhood U of 9 and constants K and M such that 

P^(e^ ^V)< Mee-^/'' 

for ip £ U. 

Proof. For n = 1 this is a simple consequence of Proposition 2, and the 
proof for general n follows by a simple inductive argument. □ 

Proof of Theorem 3. By (Dl), for any 5i > we can choose an open 
set V2 with P CV2C Uf=i^(ei,(5i) C S\D so that d{f{9), V^) > 6 for all 
9 EV2 and some constant 5 > 0. Let Vi = S\V2. Then, similarly to the proof 
of Theorem 2, we can use the decomposition 5 = Vi U V2 to write 

^11 ^12 



21 ^22 



By Proposition 2 we have 



sup P^(ef V2) < Miee-^^/'^ 
e&V2 

for some Ki and Mi, and therefore HTIiHoo ^ Miee"^^/'"^. Moreover, for 
any n > 1 we have 

(36) sup P^(e^ V2) < nMiee-^^'^\ 
e&V2 

Next we estimate ||(7Yi)'^||oo for sufficiently large A^. Using (D3), we can 
write the compact set Vi as the disjoint union U yli U • • • U A^^i U B 
where Ai = f~'^{D). (The disjointness of the Ai follows from the fact that 
f^{9) ^ D ioT 9 e D and j > 1.) Notice that for 9 e B we have r{9) D 
for all i > 0, and that E C B. 

For each 9 £ B, the condition (D2) implies there exists ng such that 
f^0{9) G V2 and f'{9) D for < z < - 1. By Lemma 6 with V = V2 
there exist a neighborhood U$ of 9 and constants Kg and Mg such that 

(37) sup F^^ie' V2) < Mgee-^'''\ 

Now suppose 9 £ Aq = D. Then f{9) and f*{9) are both in B. By Propo- 
sition 3 there exist a neighborhood Ug of 9 and constants Kg and Mg such 
that 

(38) sup P^(e^ ^ ([/;(,) U Uj.^g))) < Mgee-"^"/'". 
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Define ng = max(nj(5i), nj.(5i)) + 1. Combining (38) and (37) and (36) gives 
(39) sup P'^(e^^ V2) < Meee-^o/^^ 



for some Kg and Mg. Finally suppose that 6 £ Ai for i > 1. Then f^{9) D 
for < j < i and f^{0) S D. By Lemma 6 with V = Up(^g^ there exist a 
neighborhood Ug of 6 and constants Kg and A/51 such that 

(40) sup P*(G^ Uf^g)) < Mgee~^o/'\ 
teUg 

Define ng = i + ^p^^gy Combining (40) and (39) gives 

(41) sup F'^iQ' V2) < Mgee-^o/'\ 

for some Kg and Mg. 

We have constructed an open cover {Ug : 9 £ B U AqU • • • U A^^i} of the 
compact set Vi. Passing to a finite subcover {Ug^} and letting N = maxj ng. 
and using (36) with n = N — ng. together with (37) and (39) and (41) gives 

(42) sup P*(G^ ^V2)< Mee-^/^^ 

for some K and M, and thus ||(rfi)^||oo < Mee'^/^^ . 

It remains only to describe the eigenvalues of T22 and this can be done 
using exactly the same methods used for Theorem 2 and Remark 3. □ 
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Fig 2. Examples of the behavior described in (C). In both cases 'y — 1, I{t) = 1.4, 
g{t) = 1 + _Bsin27rt and h{t) = 0. Case (i) on the left has B — .0629, to = —.2527 giving 
f ho) = f*{to) = 1.0251. Case (ii) on the right has B = .3, to = -0863 giving f (to) = .8087, 
f*ito) = 1.473. 
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Fig 3. Plot of f for I ^ l,k = 0.1,7 = 1/12.8 
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Fig 4. Plot of f for I = l,k = 0.35, 7 = 1/12.8 
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0.2 0.4 0.6 0.8 1 
Fig 5. Plot of f (left) and p (right) for 1 = 2, k = 0.2, 7 = 1/12.8 
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0.2 0.4 0.6 0.8 1 
Fig 6. Plot of f (left) and p (right) for 1 = 2, k = 0.5, 7 = 1/12.8 
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0.2 0.4 0.6 0.8 1 
Fig 7. Plot of f (left) and third iterate f^ (right) for I = 2,k = 0.8, 7 = 1/12.8 
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Fig 8. Plot of f (left) and fourth iterate (right) for I = 2,k = 0.9, 7 = 1/12.8 
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